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We present results on the light hadron spectrum and quark mass in two-flavor QCD calculated with small sea 
quark masses down to raps/mv = 0.35. The configurations are generated using the RG improved gauge and 
tadpole-improved clover quark action a± f3 — 1.8, where ~ 1 GeV. Chiral extrapolations are made using not 
only polynomials and ChPT in the continuum but also formulae of Wilson chiral perturbation theory (WChPT) 
including 0{a^) chiral breaking terms. We examine the viability of WChPT and its influence on quark masses. 



1. Introduction 

Systematic errors associated vifith chiral ex- 
trapolations present a serious obstacle to pre- 
cise QCD calculations ^J. A rapid increase of 
the computational cost toward light quarks lim- 
its quark masses to the values larger than their 
physical ones, especially for dynamical Wilson- 
type quarks. This limitation causes sizable ambi- 
guities in the extrapolation to the physical point. 

In this report, we analyze issues related to 
small quark masses and finite lattice spacings in 
the chiral extrapolation using several chiral ex- 
trapolation functions. Our study is based on 
hadron mass data for small sea quark masses in 
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the range nips/mv = 0.60 — 0.35 |2, combined 
with our previous data at heavier quark masses 
at mps/mv — 0.80 — 0.55 jH]. These hadron mass 
data are generated with an RG-improved gauge 
action and a meanfield-improved clover quark ac- 
tion at /3 = 1.8 (a ~ 0.2 fm). We use 12^ x 24 and 
16^ X 24 lattices with the spatial size of 2.4 and 
3.2 fm. Our simulation parameters are summa- 
rized in Tabled For further details, see Ref. 0. 

2. Polynomial extrapolation 

The conventional polynomial extrapolation has 
a problem that polynomials are not consistent 
with the logarithmic singularity expected in the 
chiral limit. It is imperative to estimate the sys- 
tematic errors due to higher order contributions 
when the data are extrapolated using a low-order 
polynomial. Figure ^ shows that our new data 
at small sea quark masses deviate systematically 
from the quadratic fit of old data. The deviation 
is significant; applying higher order polynomial 
extrapolations to the whole data, Kc differs by 
lOcr, and the lattice spacing by 6.5ct(7%). 
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Table 1 

Simulation parameters at P 
runs are from Ref. 
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Figure 1. Chiral extrapolation of pseudoscalar 
meson mass. Open symbols show the results ob- 
tained in our previous calculation 

3. Chiral extrapolations based on ChPT 

A choice for controlled chiral extrapolations 
is to incorporate chiral perturbation theory 
(ChPT). However, the present lattice data are not 
quite consistent with ChPT in the continuum 
which predict the following quark mass depen- 
dence to one-loop order: 
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where Bq, f and A3 are parameters to be obtained 
by fits. Fitting the data over the whole range 
iTips/mv = 0.80-0.35 yields a large x^/dof ^ 




Figure 2. Simultaneous continuum ChPT fit to 
pseudoscalar meson mass and decay constant. 

70. We have to restrict the fitting interval to 
fnps I'^v = 0.60-0.35 to obtain a reasonable fit 
with /dof — 1.9, as represented in Fig.|21 This 
may be an indication that the chiral logarithm is 
visible only at mps/mv <0.40. 

Another possibility that we do not clearly ob- 
serve chiral logarithm in our data is the suppres- 
sion by explicit chiral symmetry breaking of the 
Wilson quark action; modifications of ChPT due 
to finite lattice spacings may be needed for anal- 
yses of data obtained on coarse lattices. 

Recently studies were made to adapt ChPT to 
the Wilson-type fermion at finite lattice spacings 
(WChPT) [51. Care has to be exercised in the 
order counting since formulae for observables are 
different. We employ the WChPT formulae which 
treat the O(a^) effects at the leading order 
because the existence of parity-broken phase and 
vanishing of pion mass depend on them in a crit- 
ical way. 

There are two features in these formulae worth 
emphasizing: (i) the coefficients of chiral loga- 
rithm terms receive 0{a) contributions, modify- 
ing the chiral behavior at finite lattice spacings, 
(ii) terms of the form ^ogrUquark appear, which 



are more singular than the 



rk log mq^ 



terms toward the chiral limit at a finite lat- 
tice spacing. The a? log mquark terms must be 
resummcd. The resummed WChPT formulae 
(RWChPT) proposed by Aoki read. 
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Figure 3. The resummed WChPT fit to pseu- 
doscalar meson mass and AWI quark mass. 
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where the fit parameters are Kc in Tn^^^^f., A, 

Wo, ^P, t^f^^, Aq, A3 and As^awi- As we 
show in Fig. |3| the one-loop RWChPT formu- 
lae is capable of fitting our data over the whole 
range mps/mv — 0.80-0.35. The fit result 
shows that the leading and the one-loop contri- 
butions in RWChPT converges well and the pa- 
rameter values are comparable with phenomeno- 
logical estimates. It is also found that if we 
make a RWChPT extrapolation using our previ- 
ous data at mps/mv = 0.80-0.55, the new data 
at mps/mv = 0.60-0.35 he well on the extrapo- 
lation curve. 

Encouraged by the features above of the 
WChPT fit, we reanalyzed the previous data 
at mps/mv — 0.80-0.55 on finer lattices with 
a = 0.16(/3 = 1.95) and 0.11 fm (/3 = 2.1) using 
WChPT and performed the continuum extrapola- 
tion. For the quark mass in the continuum limit, 
we find rri^in = 2 GeV) = 3.11(17) [MeV], 
which is smaller than the previous result ^ based 
on the quadratic chiral extrapolation by approx- 
imately 10%. 



4. Conclusions 

We have extended the study of the light 
hadron spectrum and the quark mass in two- 
flavor QCD|S] to smaller sea quark masses corre- 
sponding to mpg/my = 0.60-0.35. The new data 
clearly show systematic deviations from the previ- 
ous quadratic chiral extrapolation using the data 
in the range mps/mv — 0.80-0.55. Whereas fits 
with either polynomial or continuum chiral per- 
turbation theory (ChPT) fails, the Wilson ChPT 
(WChPT) that includes effects associated with 
explicit chiral symmetry breaking successfully fits 
the whole data. In particular, WChPT cor- 
rectly predicts the light quark mass behavior from 
data at medium heavy quark masses such as 
mps/rav > 0.5. 

While our initial WChPT analyses gave encour- 
aging results, further studies at smaller lattice 
spacings are needed to establish that the WChPT 
parameters show the expected lattice spacing de- 
pendence. 
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Grants-in-Aid of the Ministry of Education 
(Nos. 12304011, 12740133, 13135204, 13640259, 
13640260, 14046202, 14740173, 15204015, 
15540251, 16028201). 
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